Perfect State Transfer, Effective Gates and Entanglement Generation in Engineered 
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We show how to achieve perfect quantum state transfer and construct effective two-qubit gates 
between distant sites in engineered bosonic and fermionic networks. The Hamiltonian for the system 
can be determined by choosing an eigenvalue spectrum satisfying a certain condition, which is shown 
to be both sufficient and necessary in mirror-symmetrical networks. The natures of the effective two- 
qubit gates depend on the exchange symmetry for fermions and bosons. For fermionic networks, 
the gates are entangling (and thus universal for quantum computation). For bosonic networks, 
though the gates are not entangling, they allow two-way simultaneous communications. Protocols 
of entanglement generation in both bosonic and fermionic engineered networks are discussed. 
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I. INTRODUCTION 

Quantum state transfer, two-qubit gates, and entan- 
glement are essential in quantum information theory and 
quantum computation Jl| . Recently, there have been 
many proposals 0, H H/H B II II U OS El ex- 
ploiting the free evolution of spin networks for accom- 
plishing these tasks. The main idea is to minimize the 
spatial and dynamical control, which is experimentally 
challenging, on the interactions between qubits. Imper- 
fect state transfer over homogeneous spin chains has been 
studied 0, El for Heisenberg and XY Hamiltonians. A 
measurement-based state transfer scheme |sj has been 
suggested for dual-spin channels. Perfect state transfer 
@, state inversion pj, and graph state generation have 
been proposed for engineered spin chains in which the 
couplings between qubits are tunable. Quantum compu- 
tation using permanently coupled spin chains has been 
proposed [a S3- Furthermore, other dynamical proper- 
ties ^3 °f s P m chains and state transfer schemes ^2 
have also been studied. 

In this paper, we generalize the results in 0, 0, |(| for 
engineered networks. In Refs. 3 an d [1 two types of 
engineered networks which accomplish perfect quantum 
state transfer have been presented. These networks de- 
pend on the known properties of special functions, and 
hence the choice of the eigenvalue spectrum and the net- 
work couplings is limited. One of the aims of this paper 
is to show how one could "design" such engineered net- 
works without reference to any special functions. One 
simply has to choose an eigenvalue spectrum from an in- 
finite set of possibilities satisfying a certain condition [cf. 
Eq. iPHl ], which is both sufficient and necessary. The 
network couplings can then be found by solving a struc- 



tured inverse eigenvalue problem. As a consequence of 
this approach we note that even a single infinitely deep 
square well or a single harmonic well enables perfect state 
transfer from across a distance. Inspired by the recent 
rapid experimental development in optical lattices (see, 
e.g., an d references therein), our discussion will be 
presented in terms of fermionic and bosonic networks 
with the presence or absence of a boson or fermion at 
a site representing the and 1 states of a qubit. Certain 
spin networks are classified to be fermionic, as we shall 
discuss. 

Another aim of this paper is to show that effective two- 
qubit gates over remote qubit pairs can be constructed 
in those engineered networks. For fermionic ones, includ- 
ing spin chains, the effective gates are entangling and 
hence universal for quantum computation. For bosonic 
networks, the gates are not entangling, but they allow 
two-way communication for different pairs of sites simul- 
taneously without mutual interference. Finally, protocols 
for entanglement generation and transfer will also be dis- 
cussed. In contrast to the scheme proposed in ^3 > these 
protocols require minimal spatial and temporal control 
on individual qubits. 



II. ENGINEERED NETWORKS 

We start with a system consisting of spinless fermions 
(or bosons) hopping freely in a network of TV lattice sites. 
In fact, the particles need not literally be spinless, but 
they all need to be polarized in the same spin state, and 
there should not be any interactions involving spin. The 
Hamiltonian is therefore of the following form: 
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where denotes nearest-neighbor coupling, u>ij is the 
time-independent coupling constant between the site i 
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and site j, and Xj represents the strength of the exter- 
nal static potential at site j. The annihilation opera- 
tors cij obey the standard (anti)commutation relations 
for bosons (fermions) and rij = aj<2j is the number op- 
erator. This model may be considered as the strong 
tunneling limit of the Hubbard model |l5j| for fermions 
and Bose-Hubbard model |l6( for bosons. In particular, 
for one-dimensional fermionic chains, this model can be 
mapped to spin chains in which spins are coupled through 
the XY Hamiltonian 



The final state becomes 



1 N 



N-l 

3=1 3=1 

(2) 

by the Jordan- Wigner transformation [Tg. Therefore, 
such spin chains will be classified as fermionic, even 
though the individual spins are distinguishable. 

Since the Hamiltonian H commutes with the total 

rij or the total z-spin op- 



number operator n to t = X) 



erator S* ot = Ylj=i a ji tne Hilbert space can be decom- 
posed into subspaces consisting of the eigenstates of ntot 
or St ot - Furthermore, as the particles are noninteract- 
ing, the eigenstates in the n-particle subspace are the 
antisymmctrized (symmetrized) products of the single- 
particle eigenstates for fermions (bosons). 



A. Quantum State Transfer 

Quantum state transfer over a network is similar to 
the quantum random walk problem, where a variety of 
networks are equivalent to one-dimensional chains |3lll8j. 
Therefore, we will now focus on a chain of N sites. For 
j = 1, 2, AT, let \j) be the state where a single fermion 
(or boson) is at the site j but is in the empty state |0) for 
all other sites and |0) be the vacuum state where all sites 
are empty. For spin chains, |0) corresponds to the state 
where all the spins are in the spin-down state |J.) and \j) 
corresponds to a spin- up state ||) for the jth spin and 
spin-down for all other spins. The Hamiltonian in this 
single-particle subspace can be written in a tridiagonal 
form, which is real and symmetric: 
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The quantum state transfer protocol involves two 
steps: initialization and evolution. First, a quantum 
state a |0) + (3 |1) to be sent is encoded at site x. The 
initial state of the network is described by \<p x ) = 
a |0) + [3 \x). Then, the network couplings u)j and Xj are 
switched on and the whole system is allowed to evolve 
under U (t) = exp (—iHt) for a fixed time interval t = r. 
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where f^ x (r) 



-iHt 
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j\e '"' \x). Any site y is in a mixed 

fy.x( T )\ < 1) which also implies that the state 
transfer from site x to y is imperfect. Our goal here is to 
find a set of uij and Xj to realize perfect state transfer. 

In 0, it is shown that when the couplings are chosen 
such that Hn = S x , where S x is the x component of the 
spin operator of a spin S = (N— 1)/2, or alternatively 
i/jv = L • S. subject to the constraint that the z com- 
ponent of the total angular momentum J = L + S being 
zero, then a mirror inversion of eigenstates with respect 
to the center of a linear chain can be implemented. This 
implies that a quantum state at site x can be transferred 
perfectly to its mirror- conjugate site x = N—x+1, 



U(T)\<p x ) = a\Q) + e- i *»p\x) 



(•5) 



where in general 4>n ^ and a single-qubit operation on 
the site x is required to remove it, in order to reconstruct 
the original state there. 

In the next section, we will introduce a systematic way 
to find the sets of 10 j and Xj for state inversion and hence 
perfect state transfer even if we did not use any of the 
above examples. Instead of solving the eigenvalue prob- 
lem, we will first choose a desired eigenvalue spectrum [cf. 
Eq. (|12|l ] for Hn and the solutions for uij and Xj can be 
found from the spectrum and the symmetrical properties 
of Hn . It is therefore an inverse eigenvalue problem. 



B. Symmetrical properties of Hn 

State inversion by free evolution crucially depends 
on (a) the reflection symmetry and (b) the eigenvalue 
spectrum. By reflection symmetry, we mean for j = 
l,2,...,LA/2j; 



Xj - Xj 



and 



(0) 



Thus, Hn has double symmetries (also known as persym- 
metric), along both the main diagonal and the second 
diagonal. We shall now show that if the above symme- 
tries are present in Hn, then one only needs the eigen- 
value spectrum to satisfy a certain condition [cf. Eq. 
(|12fl ] in order to achieve state inversion. This condition 
will later be shown not only sufficient but also necessary 
for state transfer in mirror-symmetric networks. As a 
consequence of the above symmetries, the eigenvectors 
l e fc) = Y^j=i a j \J) nave definite parities — i.e. being ei- 
ther even or odd with respect to the mirror-conjugate 
operation j — > j. In fact, the eigenvectors can be de- 
termined [cf. Eq. H10|) ] explicitly However, we need to 
know which one would change sign when inverted. This 
can be determined by the interlacing property described 
below. 
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Let Pn(E) = Ylk=o (E—E k ) be the characteristic 
polynomials of and denote the jth leading princi- 
pal minor (i.e., the characteristic polynomial obtained by 
the first j rows and columns of a matrix) of the matrix 
EI — Hiq by Pj (E) , where I is the N x N identity matrix 
and E is a real number. With Pq = 1 and Pi = E—Xi, 
the sequence of Pj(E) is a Sturm sequence and for 
j = 1, 2, N, it satisfies a recurrence relation 



Pj(E) = (E—Xj) Pj-i(E) — w? 1 P J -_ 2 (£) 



(7) 



The Sturm sequence has an important property: the 
roots E 3 k of Pj interlace those of Pj-i — i.e. 

Ei j _ l <Ei j zl<Ei_ 2 <---<E{<Ei- 1 <Ei . (8) 



This implies that 

e®i[P N -i{E k )] = (-1) x sgnlPjv-iOEfc-i)] 



(9) 



where E k = E k and P^{E k ) — 0. We shall show im- 
mediately that this interlacing property of the Sturm se- 
quence determines the parity of the eigenvectors. 



Parity 



It is known that the coefficients a k , j = 2,3, ...,N, of 
the eigenvectors are given 0] by 



U)lU) 2 ■ ■ -OJj-l 
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(10) 



normalization condi- 



with determined by the 

tion YljLi \ a j\ = 1- We note that the par- 
ity of the eigenvectors can be determined by 
checking the relative sign of any pair of mirror- 
conjugate coefficients. For convenience, we consider 
sgn [a k N /a\] = sg-a,\P N -i(E k )/uiu 2 ■ ■■u N _ 1 ] = {-l) u x 
sgn[Pjv-i (£*)], where (-l) u = sgn [wiw 2 • • • Ww-i]. 
From Eqs. © and l |10| l. if the eigenvectors are ordered 
in decreasing eigenvalues — i.e., Eq > E\ > ■ ■ ■ > E^_i — 
the parities of them change alternatively. Since 
Pn~i{Eq) > 0, the parity of the highest energy eigen- 
state |e ) is only determined by (— 1)". It is even (i.e., 
v = 0), if all ujj > 0. As the parity changes alternatively, 
once the parity of |eo) is known, the parities of all other 
eigenvectors can be inferred immediately. These can be 
summarized as 



\e k ) = (-l) k+v \e k ) 



(11) 



for k = 0, 1, N - 1, where \e k ) = £f = i a) \j) 



C. Mirror Inversion 

Next, we require, for some time interval r, the eigen- 
value spectrum of Hjq to satisfy the relation 



-\) ±k e~ i ^ N 



where 4>m is independent of k and the ± sign has to be 
taken consistently for all k. For simplicity, we assume all 
Uj > 0. Consider U(t)\x) = E k =o e~ lEkT \e k ){e k \x). 
When \e k ) is replaced with (— l)~ k \e k ), together with 
the relation in Eq. i|12fl . the completeness relation 
/ = ^2 k= Q |efe)(efc|, and the double inversion relation 
( e k\j) = (&k\j), one can show that 



U(r) 



(13) 



From Eq. (JSJ), consequently, quantum states can be 
transported from any site x to its mirror- conjugate site 
x after a fixed period r. Once a spectrum is determined, 
the search for the solutions of ojj and Xj becomes an 
inverse eigenvalue problem. There are some efficient al- 
gorithms available in the literature for accomplishing the 
task — for example [2(j and references therein. 

Here we also note that the condition in Eq. (|12fl is 
not only sufficient but also necessary for perfect state 
transfer in mirror-symmetric networks. To prove that it 
is necessary, we set for some time r, 1 = \ (x\ U(t) \x)\ = 



Efc \( e k\x) 



-iE 



< J2k \ ( e k\ x )\ = 1, where 



(-if and we have used the normalization condi- 
tion in the last step. As the above equality must hold, 
e tlfk should be a constant phase (independent of k), and 
hence the condition in Eq. 1|12|) follows. 



Example 

Two types of spectrums, E k = —k and E k = k(k + q) 
for some rational number q and k = 0, 1, 2, N, sug- 
gested in Q can easily be shown to satisfy Eq. JT^Jl. 
However, these spectra are related to some known exam- 
ples of special functions. To illustrate the generality of 
the method, we consider a 4 x 4 tridiagonal matrix with 
the spectrum E = 1, E\ = 2, E 2 — 3, and E 3 = 2(1 +m), 
for any integer m > 1. The condition i|12|) is satisfied 
with r = 7T and = 0. One of the solutions for the 
Hamiltonian [of the form of Eq.©] is found to be 



/ a c 

c b d 

d b c 

\ c a 



(14) 



(12) 



with a = 2 + l/(2m), b = m + 2 - l/(2m), c = 
^/l — l/(4m 2 ), and d — m. The generality of generating 
engineered chains for perfect state transfer is thus clear. 
On the other hand, it is interesting to note that in the 
limit m 3> 1, one may want to put aw 2, i»mc«l, 
and d = m. However, since a ~ O(l), although m 3> 2, 
changing b from m + 2 to m would cause a large error. 
This is also confirmed numerically. Therefore, in such a 
limit, the requirement of precision is very high. In this 
sense, energy spectra that yield more uniform coupling 
are more desirable from the engineering point of view. 
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FIG. 1: In engineered fermionic and bosonic chains, when 
the spectrum of the single-particle Hamiltonian in Eq. © 
satisfies the relation 1121 . effective gates on mirror-conjugate 
pairs of sites, such as (x,x) and (y,y), can be constructed by 
free network evolution. 



Continuous Systems 

An interesting consequence of Eq. (|12[1 can also be 
found in infinite dimensional systems. The eigenvalue 
spectra allowed by Eq. I|12(l correspond to some canon- 
ical systems such as harmonic well or infinite square 
well. In those cases, the necessary criterion for state 
inversion — namely, the parity of the eigenstates — is au- 
tomatically satisfied. For example, the energy spectrum 
of an infinite square well is quadratic £7. oc fc 2 , for 
k = 1,2,3,... and the eigenfunctions ipk(x) have a def- 
inite parity ipk(x) = ( —V) k ~ 1 il)k.{x). One can show that 
(also mentioned in |21|) any single-particle wave- function 
ty(x,t) at x will be transported (up to a — sign) to — x, 
^{x,t) = (-x,t+r), for a period of r = 2-kH/Ei, 
where E\ is the ground state energy. This property has 
been discussed recently in the literature on fractional 
wave- function revivals p^ . but its relevance to quan- 
tum communication and its connection to the above gen- 
eral theory linking eigenvalue spectrum to perfect state 
transfer has not been appreciated. For example, one can 
think of the following strategy for communicating per- 
fectly through those continuous systems. We can encode 
the information of a qudit (not necessarily qubit) to the 
spin degree of freedom of a boson or fermion in a state 
^(x, t = 0) which is initially localized around x. At t = r, 
the particle will arrive x and the information can be ex- 
tracted. 



D. Effective Two-Qubit Gates 

Two identical fermions (bosons) at sites x and y is 
described by the antisymmctrizcd (symmetrized) prod- 
uct state \xy) = (1/V2) (\x)\y) ± \y)\x)), with - (+) 
sign for fermions (bosons). Similarly, the two-particle 
eigenstates \eu) = (1/V2) (|e fc )|e ; ) ± |e;)|e fc )) are also 
antisymmctrizcd or symmetrized accordingly. For many 
fermion excitations, the states are more convenient to 
be represented by a Slater determinant. Besides, in the 
mapping of spin states to fermionic states, we adopt the 
convention 0, Q that the site indices are arranged in 



ascending order. Thus, in a spin chain having spin-up 
states at x and y but spin-down for all other spins, 
the spin state corresponds to the fermionic state \xy) 
if x < y and \y x) if y < x. 

By effective gate, we mean the configuration of all the 
sites after the network evolution is the same as before, 
except that the state of the pair of qubits at x and x is 
changed according to a logic gate U x . The simplest way 
is to choose all other sites to be empty (or all other spins 
being the spin-down state for spin chains). Using similar 
tricks as before, one can show that 



U(t) \xy) = e 



— 2i<j)i\ 



\xy) 



(15) 



If y = x, then there is an extra factor (—1) for fermionic 
states but not bosonic states, after exchanging the site in- 
dices. We define a new basis with |00) = |0), 1 10) = \x), 
1 01) = \x), and |11) = \xx). In this basis, an effective 
two-qubit gate U x can be constructed readily for any con- 
jugate pair of sites x and x: 



U x = 



( 1 









? — i<t>N 










(-l)*e 



\ 



— 2i<pN 



(16) 



where 77 = 1(0) for fermions (bosons). The effective gate 
for a three-s pin chain in is a special case of U x here. 
It is known |24j that any two-qubit gate that can create 
entanglement between two qubits is universal for quan- 
tum computation, when assisted by one-qubit operations. 
Here we assume such one-qubit gates are available and we 
will show that U x can create entanglement for fermionic 
chains (and spin chains) and hence is universal for quan- 
tum computation. For bosonic chains, however, U x is 
not entangling but it allows two-way communication — 
i.e., transfer states from both ends simultaneously. 



E. Entanglement Generation and Communication 

The entanglement generation protocols in [|| can now 
be generalized. These protocols require minimal spatial 
and temporal control of the individual qubits and are 
also advantageous in that, after extracting the entan- 
gled states at sites x and x, the whole procedures can 
be repeated by replacing the extracted state with the 
corresponding initial states. The configuration of the in- 
termediate sites or spins will not be changed after each 
cycle (except the middle site/spin in protocol 1, which 
can act as a trigger of the evolution). Moreover, these 
protocols can be deployed for studying the dynamics of 
entanglement flow |23|. 



Entanglement Generation Protocol 1 

For a linear bosonic or fermionic chain with odd num- 
ber of sites, the entanglement generation problem can 
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be mapped to the state transfer problem. In 6], only 
one specific type of mapping is discussed — namely, the 
one proportional to S x . Here, with the enlarged set of 
choices for the coupling constants, we can generalize the 
mapping by including the possibility of nonzero diagonal 
coupling terms. For the sake of comparison with proto- 
col 2, we outline briefly the mapping below. 

Suppose the coupling constants still satisfy the sym- 
metry relations in Eq. |JB}, we consider a basis consist- 
ing of maximally entangled states |j) = ^ + 
for j — 1,2, ...,ti — 1, and a state \h) = |n), where 
n = h(N+l) is the position of the middle site. The 
Hamiltonian in Eq. Q acts in this basis as 



/ Ai 


U>1 








\ 




A 2 














A 3 ••• 















-1 


V o 





V2u>„_i 


A)) 


/ 



(17) 



which is also a real and symmetric tridiagonal matrix as 
Hn (but the size is about half of it). Suppose the initial 
state is |n) — i.e., a single boson or fermion at the middle 
site but empty for all other sites. The task of entangle- 
ment generation for the remote pair of sites located at 
1 and 1 is the same as to rotate from the unentangled 
state |n) to the entangled state |l). This is equivalent to 
the state transfer problem we have discussed and can be 
solved in exactly the same way. 

On the other hand, for linear chains with an even 
number of sites, a similar protocol [(| can be used for 
transferring entanglement from (now redefined) \n) = 
-^=(\n) + \n)), where n — N/2, to the remote pair of 

sites |l). However, this requires the local pair of sites 
in the middle to be maximally entangled initially and is 
therefore an entanglement transfer protocol. 

The two protocols above require the initialization to be 
made in the middle of the chains. After the free evolu- 
tion, the entangled states are then extracted at the ends 
of the chains. In situations where we are allowed to have 
access only to the pair of sites we want to entangle, proto- 
col 2, as we shall see next, will be more useful. However, 
protocol 1 works for both fermionic and bosonic chains 
but protocol 2 is applicable for fermionic chains only. 



Entanglement Generation Protocol 2 

We now show that any pair of mirror-conjugate sites x 
and x can be maximally entangled with the application of 
U x and the state initialization at x and x only. For sim- 
plicity, all other sites are set to be empty (or spin-down in 
applying to spins chains) . First of all, for any normalized 
pure state of two qubits, a |00) + b |01) + c |01) + d |11), 
where \a\ z + \bf + \cf + \d\ z = 1, the concurrence 
C = 2 \ad— bc\ is a measure of entanglement |25|. The 
two sites are unentangled when C — and maximally 



entangled when C = 1. Suppose the two sites are ini- 
tially in a product state — i.e., ad = be — and all other 
sites being empty. With the application of U x , the con- 
currence becomes 2 \ad — (—l) v bc\ = 2(1 — (— l) v ) \ad\. 
Consequently, for fermionic chains (with 77=1), the sites 
x and x can be maximally entangled from any initial 
product state with ad — be and \ad\ = j. For example, 
if the initial state is |+> |+) where |+) = (|0> + |1)), 
then the entanglement of the final state can be made ex- 
plicit by expressing it in the Schmidt form -L(|O)|0+) + 

e"^" 1 1) !</>_)), where \<j>±) = (|0) ± e~ #N |1)) and 
(</>+]</>-) =0. 

Last, we note that the protocol for generating a class of 
multipartite entangled states, called Graph states, sug- 
gested in [jj, can also be extended for the more general 
Hamiltonian in with various spectrums. 



Two-way Communication 

For bosonic chains, one can show that the entangle- 
ment of any pure state between sites x and x is invariant 
after the application of U x . In fact, the net effect of 
the free evolution of the network, with any initial con- 
figurations, for a period of t is an inversion of quan- 
tum states about the center of the chain, apart from 
an extra induced phase e~^ N . Nonetheless, this im- 
plies the possibility of simultaneous transfer of quantum 
state from site a; to site x and vice versa. Let us define 
the protocol more clearly. Suppose Alice and Bob are 
sending their states at x and x, respectively. We con- 
sider the initial state is in a product state, which can be 
written in general as (a \0) x + ai \l) x ) (b a \0) x + bi |1) 2 ), 
with all other sites being empty. Applying U x yields 
(bo \0) x + e-^h |1)J (a |0) s + e~^ ai There- 
fore, both states can be sent simultaneously. Interest- 
ingly, different parities can use the same channel, but on 
different conjugate pair of sites, at the same time without 
mutual interference. 



III. CONCLUSION 

We have demonstrated how to perform quantum state 
transfer and construct effective two-qubit gates in engi- 
neered networks in which the coupling constants are de- 
termined by the eigenvalue spectrum satisfying a certain 
condition. This condition is shown to be both sufficient 
and necessary in mirror-symmetrical networks. The pos- 
sibility of perfect communication between distant sites of 
a single harmonic trap or an infinitely deep square well 
has been discussed. The effective gates for fermionic net- 
works, including spin chains, are entangling and hence 
can be used for universal quantum computation. Two 
entanglement generation schemes are proposed. The first 
one works for both fermionic and bosonic chains but the 
second one is for the fermionic chains only. Noncthe- 
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less, the bosonic chain allows two-way communication sions and D. Burgarth for pointing out a relevant refer- 

for different pair of sites simultaneously without mutual ence. 

interference. 
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